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Abstract. We report the discovery of fast, frequency-dependent intensity variations from the scintillating intra-day variable 
quasar J1819+3845 at All cm which resemble diffractive interstellar scintillations observed in pulsars. The observations were 
taken with the Westerbork Synthesis Radio Telescope on a dozen occasions in the period between Aug 2002 and Jan 2005. 
The data were sampled at both high temporal and high frequency resolution and have an overall simultaneous frequency 
span of up to 600 MHz. In constructing the light curves and dynamic spectra the confusion from background sources has 
been eliminated. The timescale (down to 20 min) and the bandwidth (frequency decorrelation bandwidth of 160 MHz) of 
the observed variations jointly imply that the component of the source exhibiting this scintillation must possess a brightness 
temperature well in excess of the inverse Compton limit. A specific model in which both the source and scintillation pattern are 
isotropic implies a brightness temperature 0.5 x lO'^Zpc K, where previous estimates place the distance to the scattering medium 
in the range Zpc = 4-12 pc, yielding a minimum brightness temperature > 20 times the inverse Compton limit. An independent 
estimate of the screen distance using the 21 cm scintillation properties alone indicates a minimum screen distance of z a; 40 pc 
and a brightness temperature above 2 x 10'* K. There is no evidence for anisotropy in the scattering medium or source from the 
scintillation characteristics, but these estimates may be reduced by a factor comparable to the axial ratio if the source is indeed 
elongated. The observed scintillation properties of Jl 8 19+3845 at 21 cm are compared with those at 6 cm, where a significantly 
larger source size has been deduced for the bulk of the emission by Dennett-Thorpe & de Bruyn (2003). However, opacity 
effects within the source and the different angular scales probed in the regimes of weak and strong scattering complicate this 
comparison. 
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1. Introduction 

The presence of intra-day variability (IDV) in some ac- 
tive galactic nuclei (AGN) at centimetre wavelengths 
JHeeschen 19841 IWitzel et al. 1986 ) raises concern that the 
brightness temperatures of radio sources may violate the 
inverse Compton limit by several orders of magnitude 
( IQuirrenbach et al. 199'T|lKedziora-Chudczeret al. 1997> . 

It is now recognized that the variability is largely, if not 
exclusively, due to scintillation in the interstellar medium of 
our Galaxy (e.g. [Jauncey et al. 2000| iLovell et al. 2003j . This 
is established unequivocally for the IDV quasar J1819-I-3845. 
A ~90s time delay in the arrival times of the source's inten- 
sity variations measured between two widely-separated tele- 
scopes firmly identifies its variability with interstellar scintil- 
lation. The finite delay is attributed to the finite speed with 
which the scintillation pattern moves transverse to the line 
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of sight (Dennett-Thorpe & de Bruyn 2002). An annual mod- 
ulation in the timescale of the variability is also observed 
in J1819-I-3845, and is explained by the annual modulation 
in the scintillation velocity, Viss, due to the Earth's chang- 
ing orbital velocity relative to the interstellar scattering ma- 
terial ( Dennett- Thorpe & de Bruyn 2003|. This annual cycle 
arises because the Earth's velocity is comparable to the ve- 
locity of the scattering material in the interstellar medium 
responsible for the intensity fluctuations. Annual cycles are 
also reported in several other IDV sources (Bign alletal. 2003| 
IRickett et al. 2001l|Jauncey & Macquart 2001^ . 

The brightness temperatures of IDV sources have proven 
difficult to constrain, largely because of the difficulty of de- 
termining the distance to the scattering material, z, responsi- 
ble for the intensity fluctuations. These brightness tempera- 
tures have hitherto been estimated on the basis of the mea- 
surements made at A6 cm, at which flux variations are caused 
by weak interstellar scintillation. In this regime one mea- 
sures the source variability timescale relative to the Fresnel 
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timescale, fp - rp/viss - (^■z/27r)'^^/wiss and thus de- 
duces the source size relative to the angular scale subtended 
by the Fresnel scale, - rp/z- Despite these difficulties, 
several sources still exceed the inverse Compton limit by 
a large margin. A thorough analysis estimates the redshift- 
corrected brightness temperature of PKS 0405-385 at 5 x 
10'^ K ( Rickett, Kedziora-Chudczer & Jauncey 2002 1. A sim- 
ilar but less robust limit of 7"^ > 5 X 10^^ K is also derived for 
the source PKS 1519-273 ( |Macquart et al. 2000| l. 

J1819-I-385 is estimated to possess a modest brightness 
temperature of approximately ~ 10'^ K at 4.9 GHz (Dennett- 
Thorpe & de Bruyn 2003). The flux density of the source varies 
on timescales as short as 40 minutes at this frequency. These 
variations are also interpreted in terms of scintillation in the 
regime of weak scattering ( |Dennett-Thorpe & de Bruyn^ OOOi. 
The estimated screen distance is small, in the range z - 4- - 
12 pc, and the source is relatively weak, which accounts for the 
small brightness temperature despite its extremely rapid vari- 
ability f Dennett-Thorpe & de Bruyn 2003^ . 

This source generally exhibits slower variations at lower 
frequencies (jbennett-Thorpe & de Bruyn 2000l. The dominant 
variations at 2.4 and 1.4 GHz occur on ~ 6-hour timescales 
during the same interval in the annual cycle of J1819+3845 
in which dramatic intra-hour variations are observed at 
4.9 GHz. Such a change in the character of the variabil- 
ity with decreasing frequency is typical of intra-day variable 
sources JKedziora-Chudczer et al. 19971 [Macquart et al. 2000| 
IQuirrenbach et al. 20001. These slow variations are attributed 
to the increase of scattering strength with frequency, and are 
associated with refractive scintillation in the regime of strong 
scattering. Such scattering occurs when the Fresnel scale rp ex- 
ceeds the diiTractive scale length, rjiff, the transverse length 
scale over which the mean square difference in phase delay 
imposed by plasma inhomogeneities in the ISM is one radian. 
Refractive variations occur on a timescale ~r^/rdis = /"ref- 

A source scintillating in the regime of strong scattering may 
also exhibit very fast, narrowband intensity variations due to 
diffractive scintillation. These are routinely observed in pulsars 
(e.g. Rickett 1970; Ewing et al. 1970). Diffractive scintillation 
is only observable for source sizes 0,, < r^ifflz- This angu- 
lar size requirement is so stringent that no extragalactic radio 
source has previously been observed to exhibit diffractive scin- 
tillation (Dennison & Condon 1981 , Condon & Backer 1975 1. 
However, when present, it is identifiable by the fast, narrow- 
band character of its variations, which are distinct from the 
slow, broadband variations exhibited by weak and refractive 
scintillation. 

In this paper we present the discovery of narrowband, fast 
scintillation in the quasar J1819H-3845 ai All cm, characteris- 
tic of diffractive scintillation, and derive the properties of the 
source component undergoing this effect. Technical issues re- 
lated to the reduction of data from this variable source are dis- 
cussed in the following section. In Sect. 3 we derive the char- 
acteristics of the variability and derive the interstellar medium 
and source parameters associated with the phenomenon. The 
implications of the discovery are discussed in Sect. 4 and the 
conclusions are presented in Sect. 5. 



2. Observations and Reductions 

Since its discovery in 1999 J1819-I-3845 has been observed 
very regularly with the Westerbork Synthesis Radio Telescope 
(WSRT). Most observations were obtained at 4.9 GHz but since 
2001 we have increased our monitoring at 1.4 GHz as well. 
Initially 1.4 GHz observations were done with the continuum 
backend which provides 8 contiguous bands of 10 MHz cen- 
tered around 1380 MHz. As of July 2002 we have used the new 
line backend which provides both increased spectral resolution 
as well as a doubled overall bandwidth. In this paper we re- 
port on these wide-band observations that were taken on eleven 
dates between 14 Jul 2002 and 30 Jan 2005. The last three ob- 
servations were taken with an ultra wide frequency coverage 
and will be described separately below. 

2. 1 . Wideband (160 MHz) data: July 2002 till Nov 2003 

All observations were continuous over a 12-hour duration. The 
basic integration time of the WSRT is 10 s but data were av- 
eraged for either 30s or 60s. The backend was configured 
to observe simultaneously in eight 20 MHz wide sub-bands 
centered at 1450, 1428, 1410, 1392, 1370, 1350, 1330/1332 
and 1311/1310MHZ. Small gaps around 1381 and 1320 MHz 
were introduced to avoid occasional RFI at these frequencies. 
Each sub-band was further subdivided into 64 Hanning-tapered 
channels yielding 625 kHz spectral resolution. Further process- 
ing was performed only on the 28 odd channels from channel 3 
to 57. The lower and upper channels at the edges of each sub- 
band were discarded because of higher noise levels. Due to hu- 
man error no spectral tapering was applied to the 22 Feb 2003 
data. To retain full sensitivity channels 3-58 were processed for 
that epoch. The total overall bandwidth spanned by each sub- 
band is 17.5 MHz. Full polarization information was obtained 
but will be presented elsewhere. The source J1819-I-3845 shows 
only very faint (< 1 %) polarization which did not interfere with 
the total intensity analysis. 

The observations were generally taken under reasonable 
to fair weather conditions (i.e. no strong winds or precipita- 
tion). Continuous radiometry using a frontend noise source on 
a 10 second time interval provided accurate system tempera- 
tures to convert the correlation coefficients to relative flux den- 
sity. These relative flux densities are converted to absolute flux 
densities using the primary WSRT flux calibrator 3C286 which 
is tied to the Baars et al. (1977) scale, which assumes a flux 
density of 14.9 Jy for 3C286 at a frequency of 1380 MHz. The 
calibration procedure takes the spectral index of 3C286 at these 
frequencies (-0.48) into account. 

A typical example of the run of system temperatures with 
hour angle is shown in Fig. \l\ for one of the 14 telescopes 
and all 8 frequency bands. A slight non-linearity in the back- 
end leads to a small underestimation of the flux density at the 
extreme hour angles where the system temperature increases 
slightly due to increased spillover from ground radiation. We 
estimate this non-linearity effect to be about 1 % at most. (The 
non-linearity was cured in the late spring of 2004). Flux den- 
sity errors due to telescope pointing and telescope gain errors 
are well below 1% at 1.4 GHz. Overall we therefore believe 
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Fig. 1. System temperature for one telescope as a function of 
hour angle for each of the eight 20 MHz bands. Band 1 is at 
1450 MHz, band 8 at 1310 MHz. Each curve was normalized 
to an average Tsys of about 30 K. Different bands are displaced 
by 5 K. Note the effects of the radar RFI in bands 7 and 8 (see 
Section 2). 

the flux density scale to be good to 1-2%. This is corroborated 
by the relative flux density stability of the pairs of calibrator 
sources observed before and after the 12h run on J1819-I-3845 
(3C286/CTD93 before and 3C48/3C147 after). We also have 
long track observations on several known stable sources which 
agree with our 1% long-term stability assessment. 

The determination of weak, rapid intensity fluctuations in a 
radio source is not trivial when using aperture synthesis tech- 
niques at low frequencies, especially in an E-W synthesis ar- 
ray like the WSRT where 12 hours is needed to synthesize a 
good beam. It would seem to violate the principle of 'synthe- 
sis' which requires a nonvariable sky. However, there are no 
fundamental limitations in taking care of source variability, as 
is described by Dennett-Thorpe & de Bruyn (2000, 2003) for 
data taken on J1819+3845 at a frequency of 4.9 GHz. At 1.4 
GHz the situation is more complex. Within the field of view 
of the array's 25 m dishes several hundred other sources are 
detected, with flux densities, in addition to the ~ lOOmJy of 
Jl 8 19-H3845, ranging from 13mJy to O.lmJy. One of the 12h 
images is shown in Fig.|3l The noise level in the final images is 
typically 10-15 fiiy per beam. 

We have now observed this field about a dozen times at 
1.4 GHz. Although the character and magnitude of the varia- 
tions of J1819-I-3845 changes significantly, the confusion from 
the sky is expected, and observed, to be very stable, allowing 
it to be modelled well. Before vector averaging the visibilities 
for either 30 s or 60 s time intervals to form the light curve, we 
removed the response from typically 250 background sources. 
Any residual confusion from fainter sources is estimated to be 



less than 1-2 mJy which is typically 1% of the flux density of 
J1819-I-3845. More importantly, for the present study, is that 
these residual effects are broadband in nature and would be 
very similar from epoch to epoch (because the uv-coverage is 
very similar). We are therefore certain that the observed fast 
and spectral variations are not due to background confusion 
and must be due to the properties of the source and the inter- 
stellar medium. 

With observations spread across all seasons we have of 
course frequently observed in daytime. The quiet Sun still con- 
tributes a strong signal at 21 cm despite the > 40 dB attenua- 
tion by the primary beam. However, the visibility function of 
the quiet Sun drops very fast and is undetectable at projected 
baselines beyond a few 100 m. In a few cases rapid l-2mJy 
fringing was observed and short (< 144 m) baseline visibilities 
were excluded from the visibility averaging. 

The lowest sub-band is intermittently affected by interfer- 
ence due to 1 MHz of spectral overlap with a nearby radar. The 
radar, which is activated with a period of 9.6 s, beats with a 
4-minute period when observed with the WSRT, whose noise 
sources are monitored every 10 s. Close inspection of the data 
suggests that the band centred on 13 II/ 13 10 MHz is most af- 
fected, with weak interference also present on the band cen- 
tred at 1330 MHz (see e.g. Fig. El the 22 Feb 2003 dynamic 
spectrum). In the calculation of the characteristics of the scin- 
tillation signal, data from the entire two lower sub-bands are 
excluded whenever interference is evident. 

Most of the RFI in the two low frequency bands enters 
the light curves via the system temperature correction proce- 
dure. We have therefore also processed in parallel the data 
without applying this correction. The dynamic spectra for the 
1310/1311 MHz band indeed then look much cleaner To take 
care of the slight systematic system temperature variation with 
hour angle the data for February and April 2003 shown in Fig.|4] 
were corrected using the uncontaminated Tsys curves for the 
higher frequency bands . 

In order to ascertain the accuracy of the overall amplitude 
calibration on a range of timescales, we have also reduced a 
10 hour observation of the bright stable radio source CTD93, 
observed in May 2003 with a similar instrumental setup as 
J1819-I-3845. The power spectrum of the temporal fluctuations 
of this source is shown in Fig.|2l after scaling the intensity to 
a mean flux of 100 mJy. This means that the thermal noise 
has been reduced to an insignificant level and we are left with 
the combined variations due to atmospheric opacity, pointing 
and amplitude calibration. The level of 'variability' observed 
in CTD93 on frequencies of < 0.001 rad/s, which correspond 
to timescales of about 10m to 2h, is significantly less than 1%. 
On faster timescales this drops to about 0.1%, a level probably 
set by the total power amplitude calibration. All fluctuations 
observed in J1819+3845 appear to be significantly in excess of 
these levels, at any temporal scale, and become thermal noise 
hmited at the fastest timescales sampled. 
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Fig. 2. The power spectrum of temporal variations from a 10- 
hour observation of the stable bright radio source CTD 93. In 
order to compare with the temporal variations in J1819+3845, 
the power spectrum is computed from a light curve in which 
the flux densities have been reduced (by a factor ~50) to have 
a mean of 100 mJy (see §2). 



2.2. Ultra wideband (-600 MHz) data: Jan 2004 to Jan 
2005 

As of Jan 2004 the WSRT online software allowed frequency 
switching with high efficiency between different frequency set- 
tings within the L-band receiver which covers the range from 
1150 - 1800 MHz. Although a significant fraction of this 
band is aff'ected by man-made interference - GPS, GLONASS 
and geostationary satellites - the advantage of the wider band 
to study the spectral decoiTelation effects of the scintillations 
more than outweighed this loss of data. Three observations 
were taken with this ultra wideband setup, which had the 
following switching scheme: at intervals of 60 seconds we 
switched between different frequency 'combs' of 8 adjacent 
20 MHz bands. For every comb the first 20 s of data had to 
be discarded leaving 40 s of good data. The observations of 25 
Jan 2004 and 12 April 2004 were carried out with three fre- 
quency combs but diff'erent central frequencies. In the most re- 
cent data of 30 Jan 2005 we used 4 combs almost completely 
covering the available L-band frequency range. After calibra- 
tion and RFI editing the data were averaged over 40 s timeslots 
leading to Ught curves sampled on a regular 180 s or 240 s grid. 

The strong RFI encountered in several bands of each comb 
made it impossible to provide a reliable intensity calibration. 
The decomposition of the total power radiometric data into 
system temperatures and electronic gains requires stable con- 
ditions during a 10 s period, which is obviously not the case 
under strong and impulsive RFI conditions. The effects of this 
on the amplitude stability were exacerbated by the small lin- 
earity problem in the receiver. The ultra-wideband data were 
therefore internally calibrated on a band-by-band basis for each 
band of the frequency combs by normalizing on the 12-h aver- 
aged flux of J1819-H3845 itself. 
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Fig. 3. The field surrounding J1819+3845 as observed on 22 
August 2003 at 21cm. 



3. Results and Analysis 

The dynamic spectra displayed in Fig. |3 and Fig. [S] present 
a concise summary of the intensity fluctuations exhibited by 
J1819+3845 over all eleven epochs of our observations. The 
variations exhibit fine structure in both time and frequency. 
The spectral features are stochastic in nature, as both fading 
and brightening streaks are visible in all dynamic spectra. The 
spectral structure is associated with the fastest variations visi- 
ble during each epoch. This is particularly apparent during the 
22 Feb 2003 and 12 Apr 2003 observations, in which variations 
occur on timescales as short as 20 minutes, but are as long as 
several hours at other times of the year, with a variation being 
defined here as a complete oscillation in the light curve. The 
reduced duty cycle (40 s of data for every 180 s or 240 s) in the 
frequency-mosaiced dynamic spectra in Fig.|5]means that some 
of the fine temporal structure evident in other observations (c.f. 
22 Feb 2003, 12 Apr 2003) would not be as easily detectable in 
these observations. A very recent regular 160 MHz observation 
(taken on 28 Mar 2005), not presented here shows that such fine 
structure is still present. A more detailed analysis of variations 
in the 21 cm band over a period of 6 years will be presented in 
de Bruyn et al. (in preparation). 

The light curves shown in Fig. |6l further indicate that the 
intensity variations occur on several timescales and that these 
timescales change as a function of observing epoch. This is also 
demonstrated by the power spectra of the intensity fluctuations 
shown in Fig.0 The light curves from which these power spec- 
tra were computed contained no gaps over the 12 hour duration 
of the observation, with flux densities sampled every 60 s. 

We discuss the temporal variability in Sect. 13.11 and the 
spectral characteristics in Sect. l3.2l These are used to derive 
parameters of the source and scattering medium in Sect. l3.3l 
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3.1. Timescales 

Here we compare the observed power spectrum of temporal 
intensity variations to models for refractive scintillation. We 
argue that refractive scintillation fails to account for much of 
the variation observed on timescales shorter than 6 hours. 

The scintillation velocity is fastest in the period December- 
April (Dennett-Thorpe & de Bruyn 2003), so the 22 Feb 2003 
and 12 Apr 2003 datasets, which exhibit the fastest intensity 
variations, are the most useful in understanding the tempo- 
ral characteristics of the variations. Fig. 0] and Fig. |6l clearly 
demonstrate that, for both datasets, narrowband (~ 160 MHz), 
~ 20 - 120-minute variations are superposed on slower ~ 6- 
hourly variations. (The frequency-mosaiced observations made 
on similar dates in 2004 and 2005 are less suitable because 
their time-sampling is irregular, and their decreased S/N (per 
frequency channel) renders them unsuitable for characterising 
any fast, low-amplitude intensity fluctuations.) 

The variations on timescales ^ 6 hours match those ex- 
pected from refractive scintillation on the basis of observa- 
tions at higher frequencies. The transition between and strong 
scattering is thought to occur in the range ^3.5 - 5 GHz for 
this source l |Dennett-Thorpe & de Bruyn 2 00() ) and the inten- 
sity variations observed at 4.9 GHz are attributed to a scatter- 
ing screen ~ 15pc from Earth, with wiss ~ 50kms On this 
basis and assuming Kolmogorov turbulence one predicts the re- 
fractive scintillation timescale at 1.4 GHz to be between 4 and 
8 hours, consistent with the slow variations observed here. 

The larger duration of the frequency-dependent scintles ob- 
served during other epochs, their lack of associated large broad- 
band flux density deviations and the increased timescale of the 
fine, frequency-dependent structure suggests that the refractive 
scintillation timescale exceeds the 12-hour span of these ob- 
servations. This is expected on the basis of the slow-down ob- 
served in the variations at 6 cm during this period. 

Fig-Eldisplays the variations seen at 1 .4 GHz with those ob- 
served at 4.9 GHz during the same period. This demonstrates 
an excess of variability on short (20-120 minute) timescales at 
1.4 GHz. It is difficult to account for this excess in terms of 
refractive scintillation alone, since the timescale of the varia- 
tions increases sharply ( oc A^-^ (e.g. Narayan 1992; Armstrong, 
Rickett & Spangler 1995)) in the regime of strong scattering. 

To further illustrate the difficulty of accounting for all of 
the power observed in the short timescale variations in terms 
of refractive scintillation, we consider two quantitative models 
for the power spectrum of refractive variability: scintillation (i) 
from a thin phase-changing screen and (ii) from an extended 
medium in which the source size exceeds the refractive scale 

(i.e. flsi-c > 01-ef). 

The temporal power spectrum due to refractive scin- 
tillation caused by a thin screen of scattering material is 
JCodona & Frehlich 1987> 



<£/(«) = '-^ dK,(t>f^l ,kA V\ -,— 

■^iss J \viss I \vissk k j. 



X sm 



2k 



exp 



OJZ KyZ 

vissk' k 



(1) 



where o) is the angular frequency, and Ky are spatial 
wavenumbers, <1)ai^(/Ci, /c,.) is the power spectrum of electron 
density fluctuations, z is the distance to the scattering medium, 
AL <s; z is the screen thickness, is the classical electron ra- 
dius, V(r) is the visibility of the source and D^(r) is the phase 
structure function. Only the source visibility can counteract 
the sharp decline in the power spectrum due to the exponen- 
tial function at w > wiss/'"ref, but this requires a visibility that 
rises nearly exponentially quickly to account for the observed 
shallow of the decline of power spectrum (see Fig.Q. 

The refractive power spectrum for a scattering medium dis- 
tributed along the line of sight declines more slowly at high 
temporal frequencies. For a source of angular size 0s,c ^ 0ref 
scattered in a medium with thickness L and distributed ac- 
cording to Cjfiz) - C^(0) exp(-z^/L^) (with z measured from 
the observer) the power spectrum of intensity fluctuations is 
(IColesetal. 1987> 



/2 

src 



1 - exp 



\v\s,s - 



4[ 1 +(^2 ei 12] 



(2) 



For a Kolmogorov spectrum of turbulent fluctuations this 
power spectrum declines asymptotically as ~ {vi^^j ajf^^ . 

Fig-El illustrates the excess of power observed on short 
timescales relative to that expected from the two refractive 
scintillation models discussed here. The fitted models assume 
a scintillation speed wiss = 50kms"' and that all of the source 
emission is compact enough to be subject to interstellar scin- 
tillation (i.e. there are no ^ 5 mas features in the source). The 
fit parameters are listed in Table 1 . The level of fluctuations on 
timescales of less than ~ 2 hours, which the dynamic spectra 
show to be associated with highly frequency dependent varia- 
tions, are difficult to explain in terms of refractive scintillation, 
suggesting that their origin is most likely diffractive in nature. 

3.2. Spectral decorrelation of the scintillation signal 

The spectral characteristics of the scintillation signal are deter- 
mined by computing the autocovariance of the intensity varia- 
tions across the observing band, Cv(Av) = {I{v + Av)/(v)) - P. 

In order to isolate the spectral decorrelation due to the 
narrowband fluctuations only, it is necessary to remove the 
appreciable variation in the mean (intrinsic) source flux den- 
sity across our observing bandwidth. This is removed by 'flat- 
fielding' the spectrum prior to autocorrelation. The dynamic 
spectrum is normalised so that the time-average flux density 
in each spectral channel is identical, and equal to the flux den- 
sity averaged across the entire dynamic spectrum. This prevents 
the mean spectral slope across the band from masquerading as 
a scintillation signal and weights intensity fluctuations across 
the entire spectrum fairly, provided that the spectrum of the 
narrowband variations resembles the overall source spectrum. 
Although this spectral match may only be approximate in prac- 
tice, the normalisation is suflicient for present purposes since 
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22 Feb 2003 


12 Apr 2003 


z= 15 pc thin screen matching lowest frequency power 
z = 4 pc thin screen matching lowest frequency power 
z = 15 pc thin screen matching second lowest frequency power 
extended medium (AL = 15 pc) 


CjjL= 1.82 X 10'^m-i''/3 
Cj,L = 1.42 X 10"m-"/3 
C^L = 3.69 X 10''^m-"/3 
Cj, = 10-3-'8 


ClL= 1.38 X 10'^m-"'3 
ClL = 1.01 X 10"m-"/3 
ClL = 4.0 X 10" m-"'3 
Cl = 10-3-25 m-2«/3 



Table 1. Parameters used in the fits to the temporal power spectra in Fig.0 The transition frequency between weak and strong 
scattering for C^L = 1.42 x 10^^ at z = 4pc is 5.45 GHz and for z = 15 pc C^L values of 1.82 x 10"' and 3.69 x 10'^ m-'^^^ 
correspond to transition frequencies of 3.89 and 2.22 GHz respectively. 



the contribution to the error caused by any small spectral index 
mismatch is also small'. The effects of spectral misestimation 
are incorporated in the error budget when considering fits to the 
autocorrelation functions in ^3.3l below. 

The short timescale of the scintillation requires that Cy be 
computed for each one-minute time slot individually. The final 
frequency autocorrelation function is the average of the 710 - 
720 functions computed separately from each time slot. 

The autocovariances for the 25 Jan 2004, 12 Apr 2004 and 
30 Jan 2005 observations are shown in Fig.|9l These datasets 
are used because they span a sufficiently large spectral range 
that they encompass the bandwidth of the scintillation struc- 
ture. The error in our sample autocorrelation function, depicted 
by the grey region in Fig.|9l incorporates the fact that our ob- 
servations sample a finite number of scintles both in time and 
frequency (see AppendixlAi. 

Can refractive scintillation alone account for the spectral 
structure? No calculation of the form of the spectral decor- 
relation due to refractive scintillation exists in the literature, 
but refractive scintillation is known to decorrelate on a band- 
width comparable to the observing frequency: Avief ~ v (e.g. 
Goodman & Narayan 1989). We can make a simple estimate 
of the importance of refractive intensity variations by follow- 
ing the simple geometric optics model described in Narayan 
(1992). The amplitude of a typical flux variation depends on 
the root-mean-square focal length of phase fluctuations on 
scales Tref in the scattering medium. Since the focal length 
is much larger than the distance to the observer, and assum- 
ing Kolmogorov turbulence, the amplitude is ~ '"diff^^/^) 
Thus the expected flux density change across a bandwidth Av 
is A5(l + Av/v)'''/^° - A5, where AS ~ 25mJy is the root- 
mean square amplitude of the refractive scintillations. Thus the 
contribution due to refractive frequency variations is at most 
2.0 mJy across a bandwidth of 200 MHz, considerably smaller 
than the variation observed. 

Larger chromatic effects are possible if a strongly refract- 
ing 'wedge' or prism of interstellar material is present along 
the line of sight. Chromatic refraction due to a wedge would be 
important if its gradient were sufficient to change the angle of 
refraction from the low to high edge of the observing band by 



' The error is no larger than the mean square of the flux density er- 
ror; the contribution cancels to first order in flux density difference 
since the autocorrelation averages out the contributions from both 
higher and lower frequencies. For instance, the additional contribu- 
tion to the autocorrelation is ~ 0. 1 mjy^ across 100 MHz for a diffrac- 
tively scintillating component of spectral index 1 .0 and a source with 
spectral index 0.3. 



an amount comparable to 6^^f. However, if the gradient is not 
perfectly aligned orthogonal to the direction of the scintillation 
velocity one expects the scintillations to be displaced in time as 
well as frequency. The wedge would have to displace the scin- 
tillation pattern by a relative distance 0.3 r^^f from one edge 
of the observing band to the other to account for the chromatic 
nature of the intensity fluctuations observed here. One would 
then expect the wedge to also displace the scintillations in time 
by S 0.3 fief as one moves across the observing band. No such 
temporal displacement is observed. Moreover, as the scintilla- 
tion velocity changes through the year, one would expect a sys- 
tematic change in the slope of the frequency-dependent scintles 
(i.e. dv/dt) with a change in the direction of the scintillation 
velocity. As no systematic change with scintillation velocity is 
observed, we conclude that there is no spectral contamination 
due to a refracting wedge. 

3.3. Diffractive scintillation source characteristics 

The temporal and spectral characteristics of the intensity vari- 
ations are combined to determine the parameters of the source 
undergoing diffractive scintillation. We concentrate on fits to 
the spectral decorrelation of the 25 Jan 2004, 12 Apr 2004 and 
30 Jan 2005 datasets where the spectral coverage exceeds the 
typical bandwidth of the scintillation structure. 

A detailed interpretation of the scintillation parameters de- 
pends in detail on whether the scattering material is located in 
a thin layer or whether it is extended along the ray path. The 
manner in which the scale size of the scintillation pattern and 
the form of the spectral decorrelation are altered by source size 
effects depends on the distribution of scattering material along 
the line of sight. We consider two specific models, one in which 
the scattering material is confined to a thin screen a distance z 
from the observer, and the other in which the scattering mate- 
rial is extended homogeneously out to a distance Az from the 
observer 

In the thin screen model the spectral autocorrelation takes 
the form JChashei & Shishov 19761 IGwinn et al.l998t 



'thin 



= Aotf -I- 5 



diff 



1 + ^ 



Ay\, 



(3) 



where 5 diff is the flux density of the source component exhibit- 



ing diffractive scintillation, Avt 



'diff 



/ is the decorrelation 



bandwidth that a point source would possess in this scattering 
medium, and ^o/^cr is the ratio of the source angular radius to 
a critical angular scale of the scintillation pattern 0cr - rml'^Z- 
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Equation Q was first derived for gaussian turbulence but it 
has also been derived approximately for Kolmogorov turbu- 
lence dGwinn et al. 1998T l. The spectral autocorrelation in the 
thin screen model is degenerate to the combination S , Avt and 
Oq/Oct, so it is only possible to fit for two of these three param- 
eters. The uncertainty in the base level of the observed spectral 
autocorrelation function necessitates the introduction of the ad- 
ditional constant Aoff. When the source size exceeds the critical 
angular scale Oa- the form of the spectral decorrelation simpli- 
fies to 

\2 



diff 



(£). 



Av ^ Avt^ 



Ay ^ Avtf^ 



(4) 



The characteristic scale of the frequency pattern is thus set by 
the source size to 



A A 

Av = Avt — 



(5) 



The scale of the observed diffractive scintillation pattern is 



(6) 



This scale can be equated directly to the product of the diffrac- 
tive scintillation timescale and velocity, viss fdiff ■ A notable fea- 
ture of the thin-screen model is that the pattern scale grows ar- 
bitrarily large with source size. 

The spectral decorrelation associated with the extended 
medium model is cChashei & Shishov 1976 ). 

Aoff -1-5 



where /(Av) 




Both models assume that 6q/0ci- is a constant and do not ac- 
count for any possible variation of source size relative to the 
critical angle, 6„, with frequency. Both also assume a circu- 
larly symmetric source and an isotropic scattering medium, the 
implications of which are discussed further below. 

3.3.1. Simple Brightness Temperature Estimate 

The size of the feature associated with the narrowband inten- 
sity fluctuations can be estimated in a simple manner for the 
thin screen model given the transition frequency and the dis- 
tance to the scattering screen. The transition frequency, v,, can 
be deduced from fits to the refractive power spectrum given the 
distance to the scattering medium and the scintillation speed 
(Fig. Hand Table 1). Dennett-Thorpe & de Bruyn (2003) ar- 
gue that the screen distance is in the range z = 4 - 12 pc. We 
estimate the brightness temperature for screen distances of 4 
and 15 pc assuming a scintillation speed of uiss - 50kms"', a 
value deduced on the basis of measurements at 6 cm (Dennett- 
Thorpe & de Bruyn 2003). 

The amplitude of the frequency-dependent scintillation 
is 20mJy and its decorrelation bandwidth is 170 MHz (see 
Fig-E)- The flux density associated with the scintillating com- 
ponent is Sv ~ 2O/(0o/6'cr)mJy. The ratio i0o/O„) is estimated 
directly from the ratio of the observed to point source decorre- 
lation bandwidth by employing eq. This equation is valid 
in the present case since the observed decorrelation bandwidth 
is several times larger than the point source decorrelation band- 
width for the expected value of v, ^ 4 GHz (see Table 1). The 
critical angular scale O^r is estimated directly from the scatter- 
ing screen distance and the transition frequency. One solves for 
6o using Ocr and the ratio 6o / 6^ to obtain. 



cosh(V£)-.cos(V£)' '-^'-'Wc] (i^) 



(11) 



R{Av) = 2 



cosh 



Av 

AVex 



Av 

Avpv 



(8) 



where the point-source decorrelation bandwidth is Avgxt = 
TTv^r^jjj/Az. When the source is extended the spectral decor- 
relation function is nearly degenerate to a combination of the 
free parameters, but it takes the following simple form 

3 (l)' , Av « AVex 

Fex-Aoff + 5,2;^ RjA (^f , AVex ^ Av 5 AVex (|^)^ ^9) 

^ R{Av), Av » Avex (^)' . 

Source size reduces the overall amplitude of the spectral auto- 
correlation function by a factor flg/^cr- The diffractive pattern 
scale is largely insensitive to source size, and grows to a max- 
imum of only twice the diffractive scale length as the source 
size increases: 



We use the fits to the refractive power spectrum in Fig.Q 
to estimate the transition frequency. For z - 15 pc one has 
V, = 3.89GHz and raifr = 2.19 x lO^m at 1.4GHz and the ex- 
pected decorrelation bandwidth of a point source at 1 .4 GHz is 
43.3 MHz. The corresponding numbers for a screen at z = 4pc 
are v, = 5.45 GHz, raiff = 6.39 x 10*'m and Vdc = 13.8 MHz. 
For z- 15 pc the source size \s 9o - 19 fias, while for z = 4pc 
the source size is 67 juas. 

The transition frequency cancels out of the expression for 
the source brightness temperature leaving. 



(12) 



1 + 



(10) 



where we correct the brightness temperature for the source red- 
shift zs = 0.54 and we take A - 0.21 m. For z = 15 pc the im- 
plied brightness temperature is 7 x 10'^ K while for z = 4pc it 
is2x 10'^ K. 

Two factors account for the higher brightness temperature 
estimated at 1 .4 GHz relative to that at 4.9 GHz: (i) the wave- 
length is 3.5 times larger and (ii) for screen distances z ~ 15 pc, 
the source size is estimated to be up to ~ 3 times smaller than 
the value estimated at 6 cm. 
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3.3.2. Fits to the spectral autocovariance 

For completeness we also estimate the brightness temperature 
using parameters extracted from the fits to the spectral autoco- 
variance, without imposing any external constraints (e.g. from 
the fits to the refractive power spectrum). 

Fits to the frequency autocovariance functions are shown 
in Fig.|5] Both models provide close fits to the data, with re- 
duced values less than unity. Fit parameters and confidence 
limits for both thin-screen and extended medium models are 
listed in Table 2. We note that strong RFI causes difficulties in 
the spectral normalisation of the 12 Apr 2004 dataset, which 
introduces systematic errors in the frequency autocovariance 
function and the scintillation parameters derived from it. The 
scintillation and source parameters derived from this dataset 
should be treated with caution. 

In the thin screen model we fit for the products Avt^dift and 
'^diff^o/^ci, leaving the flux density of the scintillating compo- 
nent as a free parameter. In the extended medium model we 
fit to the combination S^^^bIJOq because there is a strong de- 
generacy between between 5difF and ^o/^ci- once the source size 
exceeds flcr- 

Equations Q and (|6j or Q and (llOl l are used in conjunction 
with the fit parameters and scintillation velocity and timescale 
to derive the screen distance, source size and brightness tem- 
perature for either thin- or thick-screen models. These derived 
source parameters are listed in Table 3. 

Many of the quantities derived in Table 3 depend on 
the scintillation timescale. The scintillation timescale at each 
epoch is determined by computing the intensity autocovari- 
ance function for each 17.5 MHz band and measuring the 
point at which fliis falls to 1/e of its maximum value. The 
mean timescale at each epoch is computed by averaging the 
timescales derived for the various bands. The scintillation 
timescales are 67+3 min, 40+2 min and 87±3 for 22 Feb 2004, 
12 Apr 2004 and 30 Jan 2005 respectively. The quoted error in 
the timescale is the standard error of the mean, computed from 
the variation in timescale observed between bands. The band- 
averaged temporal autocovariance functions for these dates are 
shown in Fig.|Sl In Table 3 the scintillation timescale is ex- 
pressed as a multiple of 1 hour. 

We reject the extended medium model because its estimates 
of the scattering medium properties are unphysical. The large 
medium depth indicated by this model would place most of the 
scattering medium well outside the Galactic plane. Formally, 
the estimated depth is large because, for an extended source 
scintillating in an extended medium, the diffractive pattern 
scale S() asymptotes to twice the diffractive scale length r^ig. 
This fixes the diffractive scale length to a much larger value 
relative to the thin screen model, in which the pattern scale 
increases with source size without bound. Thus, in the ex- 
tended medium model, one requires a large medium depth for 
a given decorrelation bandwidth since the latter is proportional 
to ?"jijf/z. It should be noted that just such a misestimate is ex- 
pected to occur when the source is extended but is in fact sub- 
ject to scattering through a thin screen. 

We regard the range of the brightness temperatures derived 
at diff'erent epochs as the most faithful estimate of their true 



uncertainty. Part of the range can be attributed to the uncer- 
tainty in the scattering speed and scintillation timescale, both of 
which change between epochs. The minimum brightness tem- 
perature is uncertain by a factor of two even when reasonable 
variations in these parameters are taken into account and the 
RFI-afflicted 12 Apr 2004 dataset is excluded. The uncertainty 
may reflect other uncertainties not taken into account by the 
model, such as anisotropy in the source and scattering medium. 

We note that the nearby pulsar PSR 1 8 1 3 +40 1 3 is observed 
to exhibit a diffractive decorrelation bandwidth of ~ 10 MHz at 
1 .4 GHz (B. Stappers, private communication). If the scattering 
properties of J1819+3845 are comparable then this decorrela- 
tion bandwidth favours a value of Sdm around lOOmJy, and a 
brightness temperature toward the low end of its allowed range. 

3.3.3. Source spectral changes 

In the scintillation model above the brightness temperature de- 
pends on the free parameter 5diff. We discuss here the extent 
to which the spectrum of the component undergoing scintilla- 
tion matches the mean source spectrum, and whether this com- 
ponent makes a significant contribution to the overall intrin- 
sic source spectrum at low frequencies. The latter might be 
expected if the source is comprised of multiple components 
with distinct spectra. Evidence that this may be case comes 
from analysis of the 4.9 and 8.4 GHz light curves, in which 
distinctly different polarization and total intensity fluctuations 
imply that the source is composed of at least two bright features 
(Macquart, de Bruyn & Dennett-Thorpe 2003). 

Despite the source's complex structure, its mean spectrum 
between 1.4 and 8.4 GHz is a power law with a spectral in- 
dex of 0.8. We have also measured the spectral index of the 
source from the variation in mean flux density across the band 
from our 21 cm observations. These measured spectral indices, 
listed in Table 4, varying between 0.8 and 1.2, are, with one 
significant exception, consistent with the intrinsic spectral in- 
dex derived on the basis of long-term measurements between 
1.4 and 8.4 GHz (de Bruyn et al., in prep.). It is difficult to 
be more precise, since our 1 .4 GHz spectral measurements are 
a poor indicator of the intrinsic source spectrum when only a 
few diffractive scintles are observed, as is the case in many of 
our observations. 

The one notable exception is the spectrum measured on 22 
Feb 2003, which is wildly at variance with the mean source 
spectrum and with the other observations. This difference may 
be significant, because the average spectrum extracted from 
this observation encompasses many diffractive scintles. This 
difference may reflect the emergence of a new component in 
the source. However, no discernible deviation in the amplitude 
of diffractive scintillation is associated with this epoch. 

4. Discussion 

4.1. Robustness of the brightness temperature 
estimate 

The ^ lO''* K brightness temperature implied by the diffractive 
scintiUation properties of J1819+3845 is difficult to account 
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fit parameter 


25 Jan 2004 


12 Apr 2004" 


30 Jan 2005 


thin screen 


ofi'set A„ff (Jy^) 
bandwidth AvtSae (MHz Jy) 


-3.46 ± 0.03 X 10-* 
4.09 + 0.05 
1740 ± 10 


-2.98 ± 0.04 X 10-* 
4.68 + 0.07 
1680 ± 20 


-2.93 ± 0.02 X 10-* 
3.89 ± 0.03 
1600 ± 10 


extended 
medium 


offset A„„ (Jy^) 
component flux density and source size S^gOjjOl (Jy^) 
bandwidth Avext (MHz) 


-3.10 ± 0.03 X 10-* 
1.78 ±0.01 X 10-* 
5.91 ± 0.07 


-2.60 ± 0.04 X 10-* 
1.81 ±0.02 X 10-* 
6.7 ±0.1 


-3.06 ± 0.02 X 10-* 
2.11 ±0.01 X 10* 
5.91 ± 0.06 



Table 2. Fit parameters with formal Icr errors from the fit. f Strong RFI affected the spectral normalization of the 12 Apr 2004 
observation, which in turn affected estimation of the frequency autocovariance upon which this fit is based. 



model 


quantity 


scaling 


25 Jan 2004 


12 Apr 2004 


30 Jan 2005 


thin screen 


screen distance (pc) 
source size (//as) (Sus = 0.05 Jy) 
source size (//as) (5jiff = 0.15 Jy) 
brightness temperature (K)(5difF - 0.05 Jy) 
brightness temperature (K)(5diff = 0.15 Jy) 


'ifir^SO 


6.0 
4.4 

15 

9.1 X 10'* 
2.4 X 10'* 


5.5 
4.8 
16 

7.6 X 10'* 
1.9 X 10'* 


6.9 
3.8 

13 

12 X 10'* 
3.1 X 10'* 


extended medium 


medium thickness (kpc) 
source size (juas) 
brightness temperature (K) 


5 diff tilv^o 
diff hhr^ 50 


5.7 
3.9 
1.5 X 10'« 


7.2 
3.1 
2.4 X 10'« 


5.7 
3.6 
1.7 X 10'« 



Table 3. Source parameters derived from the best- fit parameters of the various scintillation models appUed to the 25 Jan 2004, 12 

Apr 2004 and 30 Jan 2005 observations. A range of parameters are permitted by the thin screen model because the diffractively 
scintillating flux density is unknown. A larger flux density implies a greater source size and lower brightness temperature. The 
maximum possible flux density is the intrinsic flux density of the entire source, which is approximately 150 mJy. The numbers 
in the three last columns should be multiplied by the scaling parameter to derive the correct value of each quantity. Here S diff is 
measured in Jansky, tihr is the diffractive scintillation timescale in hours and vso is the scintillation speed normahsed to 50 km . 



date 


mean flux density (mJy) 


spectral index a 


14 Jun 2002 


85 


0.79 ± 0.01 


30 Aug 2002 


80 




22 Feb 2003 


144 


0.29 ± 0.01 


12 Apr 2003 


111 


1.01 ±0.01 


19 Jun 2003 


115 


1.19 ±0.01 


22 Aug 2003 


242 


(1.21 ± 0.02) 


18 Nov 2003 


151 




25 Dec 2003 


147 


1.05 ±0.01 


22 Feb 2004 


144 (at 1.40 GHz) 


0.90 ± 0.06 


12 Apr 2004 


143 (at 1.40 GHz) 


0.67 ± 0.02 


30 Jan 2005 


167 (at 1.40 GHz) 


0.50 ± 0.03 



Table 4. The variation in mean flux density and spectral index of J1819-I-3845 with observing date. Here the spectral index a, 
defined as S oc v", is derived solely on the basis of the spectrum exhibited across the band at 1.4 GHz. Blanks and bracketed 
values indicate observations in which so few diffractive scintles are present that the mean spectrum is a poor representation of 
the mean source spectrum, and blanks indicate instances for which a power law is an unacceptable fit to the mean spectrum. The 
errors quoted in the spectral index reflect only formal errors associated with a fit to a power law. 



for using the standard interpretation of AGN radio emission in 
terms of synchrotron emission. In this section we consider the 
robustness of this estimate. 

The greatest source of error in the thin-screen model is as- 
sociated with the effect of anisotropy on the scale of the scin- 
tiUation pattern, sq, which propagates into the estimation of 
fdiff. Scattering measurements of pulsars (Mutel & Lestrade 
1990; Spangler & Cordes 1998) suggest that the maximum de- 
gree of anisotropy expected due to turbulence in the interstellar 
medium is 3 : 1. However, intensity variations in the regime of 



weak scattering at 4.9 GHz indicate that the scintillation pattern 
of J1819-H3845 has an axial ratio of 14+^^° (Dennett-Thorpe 
& de Bruyn 2003). The relative contributions of medium and 
source to this overall anisotropy are unknown at this frequency. 
At 1 .4 GHz the source structure is expected to be the primary 
agent responsible for any anisotropy in the scintiUation pattern 
because the source substantially exceeds the critical angular 
scale of the diffraction pattern. We estimate 3.2 < 60/ Oct < 37 
for 50 < 5 diff < 150 mJy (see Table 2). 
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Anisotropy in the source must couple with anisotropy in- 
trinsic to the turbulence in the scattering medium to cause an 
appreciable misestimate of the source size. This is because any 
source elongation oriented parallel to the medium anisotropy 
would not be detected by a reduction in the scintillation am- 
plitude. In this case an anisotropic ratio of f would lead to a 
misestimate of the source size by a factor of ^ along one source 
axis and a brightness temperature misestimate of a factor of ^. 
Anisotropy in the source alone is insufficient to cause a serious 
misestimate of its brightness temperature, because the source 
extension would manifest itself through a reduction in the mod- 
ulation amplitude of the scintillation; although the source is 
assumed to be circularly symmetric in our fits above, the re- 
duction of the modulation amplitude would lead us to deduce a 
source angular size somewhere between the lengths of the short 
and long axes of the source. 

We have performed an analysis of the variation in scintil- 
lation time scale at 21cm. An annual cycle in the timescale 
is clearly observed, but anisotropy in the scintillation pattern 
is not required to reproduce the timescale variations observed 
from our observations to date. Changes in the magnitude of 
the scintillation velocity, as the Earth's velocity changes with 
respect to the scattering medium's, are alone sufficient to re- 
produce the annual cycle. 

Another shortcoming of the scattering models resides in the 
assumption that the scattering occurs in the regime of asymp- 
totically strong scintillation. The scattering strength rp/raift 
is derived from the decorrelation bandwidth. The decorrela- 
tion bandwidth in the thin-screen model is at most Av, = 
78 MHz for S^ff - 50mJy, which implies a scattering strength 
X! 4.2. For Sdiff = 150 mJy the decorrelation bandwidth is 
26 MHz, implying a scattering strength of 7 . The scintillation 
is sufficiently strong to be applicable to the present situation. 
Certainly, any errors introduced by this approximation are mi- 
nor relative to those introduced by possible anisotropy in the 
scintillation pattern. 

4.2. Relationship to 6 cm source structure and 
scattering properties 

It is important to consider how the structure of the source de- 
rived here relates to that inferred on the basis of the weak scin- 
tillation exhibited by the source at 4.9 GHz. 

The source angular size derived at 1 .4 GHz for any given 
screen distance (c.f. Sect. 13.3. 1> is approximately three times 
smaller than that inferred at 4.9 GHz. This would seem sur- 
prising because one would expect that the source, whose spec- 
trum does not fall off as fast as that of a uniform synchrotron 
self-absorbed source (q'= -1-2.5) should be substantially larger at 
1.4 GHz. A straightforward comparison between the scintilla- 
tion properties at 1 .4 GHz and 4.9 GHz is complicated by these 
opacity effects. Significant parts of the source that are visible 
at 4.9 GHz, and contribute to the observed scintillations, may 
well be hidden at 1.4 GHz. On the other hand, the diffractively 
scintillating component need only comprise a small fraction of 
the total source emission at 6 cm because weak and diffractive 
scintillation are sensitive to structure on widely different angu- 



lar scales. Weak scintillation responds to all structure on an- 
gular scales ^ Op whereas diffractive scintillation produces a 
strongly identifiable response specifically to structure on much 
smaller angular scales, ~ Oa-- Thus the small structure respon- 
sible for the diffractive scintillation at 1.4 GHz may well be 
present at 4.9 GHz but the signature of its presence could be 
masked by the dramatic variations due to the rest of the source. 
It is also possible that the source has multiple components at 
both 1.4 GHz and 4.9 GHz with very different spectral indices. 
For example, a coherent emitter (see the next section) could 
well have a very steep spectrum which would contribute a neg- 
ligible fraction of the emission at 4.9 GHz. 

Recent observations suggest that this may indeed be the 
case. We may have detected evidence for the component re- 
sponsible for diffractive scintillation in the weak scattering of 
the source at 4.9 GHz. Observations from Dec 2003 to Apr 
2004 at 4.9 GHz indicate the emergence of 5-10 mJy variations 
on a timescale of < 15 min superposed on the ~ 200 mJy peak- 
to-peak, ^ 40 min variations normally observed at 4.9 GHz at 
this time of year. It is possible that the scintillation at 1 .4 GHz, 
which is more sensitive to fine structure, first detected the same 
feature which was subsequently detected at higher frequencies. 
We continue to monitor the source and will return to this ap- 
parent evolution in the future. 

The screen distance indicated by the model in Sect l3.3.2l is 
larger than the 4 - 12pc value estimated by Dennett-Thorpe 
& de Bruyn (2003) assuming isotropic turbulence to model the 
intensity fluctuations observed in the regime of weak scattering 
at 4.9 GHz. The minimum distance implied by the present thin- 
screen model is ~ 40 t^i^^.v^p^pc if Sdi« = 150 mJy. An obvious 
reason for this discrepancy is that anisotropy is not taken into 
account in the estimate of the screen distance at 4.9 GHz, and 
we have no detection of anisotropy at 1.4 GHz. 

4.3. Problems with high brightness temperature 
emission 

The high brightness temperature exhibited by a compo- 
nent of J1819H-3845 raises concerns regarding the interpre- 
tation of AGN emission in terms of incoherent synchrotron 
radiation. Inverse Compton scattering limits the brightness 
temperature of incoherent synchrotron emission to 10'^ K 
(jKellerman & Pauli ny-Toth 1969] l but equipartition arguments 
( Readhead 1994 ) suggest that the actual limit should be an 
order of magnitude below this. Bulk motion with a Doppler 
boosting factor S ^ 100 is required to reconcile the ob- 
served brightness temperature with its maximum possible rest- 
frame value. Such high bulk motions are problematic because 
they imply unacceptably high jet kinetic energies. Synchrotron 
emission is also extremely radiatively inefficient in this regime, 
and it is questionable whether T S 100 motions are compat- 
ible with the hypothesis of incoherent synchrotron radiation 
( |Begelman, Rees^fe Sikora 1994> . 

In view of the difficulties confronted by an explanation 
involving synchrotron radiation, it is appropriate to consider 
whether a coherent emission mechanism provides a more ac- 
ceptable explanation of the high brightness temperature. The 
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nature of coherent emission requires the source to be com- 
posed of a large number of independently radiating coherent 
'bunches', with individual brightness temperatures far in ex- 
cess of the Tj, ^ 10'"* K value derived here (e.g. lMelrose 199ll . 
This is because the coherence volume of any one coherent 
bunch is microscopic compared to the light-week dimensions 
of the source. Further, the short lifetime associated with any in- 
dividual coherent bunch would require emission from a large 
number of independent subsources to explain the constancy of 
the emission observed on 12-hour to 6-monthly timescales. 

Coherent emission from each bunch is expected to be 
highly polarized. The upper limit of 1 % overall source polar- 
ization at 1.4 GHz limits the polarization associated with the 
diffractive component from k 1%, for S^is - 150 mJy, to 3%, 
for S diff - 50 mJy. This suggests that either the emission is effi- 
ciently depolarized as it escapes the source or that the emission 
is intrinsically unpolarized. The latter would occur if the mag- 
netic field is highly disordered within the emission region, so 
that the polarizations of individual coherent patches would be 
diluted when averaged over the entire region. 

Another important obstacle relates to the escape of ex- 
tremely bright emission from the source region. Induced 
Compton scattering places an extremely stringent limit on the 
thermal electron density of the source: for a path length L the 
electron density must satisfy 

o-tL \5x 10'^ K I \lpc/ 1 1015 K/ ^ 

for induced Compton scattering to be unimportant. It is ar- 
gued that this density is incompatible with the high densities 
required to efficiently generate coherent emission in the first 
place (e.g. |Coppi, Blandford & Rees 1993| . This difficulty may 
be overcome by appealing to a highly anisotropic photon dis- 
tribution. However, this explanation is also problematic be- 
cause such highly beamed emission acts like a particle beam 
in exciting Langmuir waves which also scatter the radiation 
(|GedaHn & Eichler 19931 ILuo & Melr ose 1995 1. This eff^ect is 
the accepted mechanism for the occultation of several eclipsing 
radio pulsars whose radiation propagates through a relatively 
low density stellar wind. Applied in the context of AGN, this 
effect would require unreasonably low electron densities in the 
emission and ambient media to permit the escape of coherent 
emission from the source (Levinson & Blandford 1995 1. 

Begelman, Ergun & Rees (2005) have recently proposed an 
electron-cyclotron maser model for the high brightness temper- 
ature emission inferred in some IDV sources. They discuss in 
much more detail the difficulties associated with the escape of 
bright radiation from a source. They conclude that it is possi- 
ble for the high brightness radiation observed in J1819+3845 
to escape, subject to certain constraints on the location of the 
emission region. 

5. Conclusions 

We have detected the diffractive interstellar scintillation from 
the quasar J1819-I-3845 at 1.4GHz. This detection is notable 
because it constitutes the first detection of this phenomenon in 



an AGN, and it implies that a component of the source must be 
extremely compact. 

These scintillations are analysed in the context of thin- 
screen and extended-medium models for the distribution of 
interstellar scattering material. The timescale, bandwidth and 
amplitude of the variations at 21 cm imply a brightness tem- 
perature > lO'^^K. 
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Appendix A: Determination of errors in the 
spectral autocorrelation function 

The finite number of scintles measured during each observation 
must be taken into account when comparing the spectral auto- 
covariance curves shown in Fig.|9]to the theoretical models, 
which represent an average over the ensemble of all possible 
scintles. 

The deviation in the sample autocovariance function from 
the ensemble-average autocovariance can be characterized in 
terms of the the variance in the sample autocovariance function, 
C(t), given by tJenkins & Watts 19681 eq. 5.3.21) 

var[C(r)] = \ f {T - t - r){y(rf - y(r + T)y(r - t)] dr, 

(T - \T\y j-(T-T) 

(A.l) 

where T is the length of the sample dataset and 7(t) is the 
ensemble-average autocovariance. 

We wish to determine the error in the spectral autocovari- 
ance of a single time slot. The ensemble average spectral auto- 
covariance is unknown but, in order to determine the eiTor in 
Cv, we approximate its form using the best-fitting thin screen 
autocovariance function. The error for a single time channel is 
determined by numerically integrating this model in equation 
JA.U . with T set to the total bandwidth spanned during the ob- 
servation. 

The final spectral autocorrelations shown in Figure|5]repre- 
sent an averages over many diffractive scintles observed during 
a 12-hour duration. We assume that the autocovariance esti- 
mates of the individual diffractive scintles are statistically inde- 
pendent, so that the errors deduced from eq. iA.H are reduced 
by the square root of the number of diffractive scintles con- 
tributing to Cv This is estimated by comparing the observing 
duration to the diffractive timescale. 

These sampling errors are combined linearly with spectral 
misestimation errors (see ^3.2> to derive the error regions de- 
picted in Fig.|5] Sampling errors dominate at all spectral lags. 
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Fig. 4. The dynamic spectrum of J18 19+3845 over 160 MHz on eight epochs in the interval Jun 2002-Dec 2003, arranged by day 
of year, given above the panels. Each point in each dynamic spectrum represents a one minute average over a 625 kHz bandwidth. 
Bad data flagged in the analysis are represented by black points. The periodic fringes visible in the lowest frequency band on 22 
Feb 2003 and 12 Apr 2003 are caused by terrestrial interference. Note that the displayed flux density range in some epochs when 
the source is not strongly variable (e.g. 30 Aug 2002) is rather small, making the dynamic spectrum appear more noisy. 
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Fig.5. Wideband, frequency-mosaiced observations of the scintillation spectrum in J1819+3845 on 25 Jan 2004 (top), 12 Apr 
2004 (middle) and 30 Jan 2005. The sizes of the white vertical spaces between data channels are proportional to their correspond- 
ing frequency gaps. Bad data points, channels or bands are marked by black pixels. The contribution of the intrinsic spectrum 
of J 1819+3845 has been removed by normalising each spectral channel so that they are identical (and equal to the mean source 
flux density, averaged over the entire dynamic spectrum). 
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Fig. 6. Light curves of the variabiUty on 22 Feb 2003 and 12 Apr 2003 at 1.4 GHz in each band. The (colour-coded) light curves 
from individual bands differ on timescales ranging from 20 mins up to several hours. The variations at the lowest frequency 
band in the 12 Apr 2003 dataset were contaminated by terrestrial interference, and are excluded from this plot. The observing 
frequency increases with colour from red to blue. The differences between the band light curves are well above the 1 mJy noise 
level (per 60 s). The bottom panel, a light curve of the source at 4.8 GHz on 2 Feb 2003, contrasts the smoothness of the intensity 
variations observed in the regime of weak scattering with those observed at 21cm. The noise on the 4.8 GHz light curve is 
typically 1 % of the flux density or less. These errors are discussed in Dennett-Thorpe & de Bruyn (2003) in detail. 
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Fig. 7. The average power spectra of the temporal variations on 22 Feb 2003 (top) and 12 Apr 2003 (bottom) shown with the 
expected thin screen and extended medium power spectra for comparison. The lines are labelled according to the model: (1) 
best-fitting extended medium with medium scale length 15pc, (2) thin screen model with z - 15 pc with C^AL adjusted to 
match the power at the lowest frequency bin (3) same as (2) but with C^AL adjusted to match the power at the second-lowest 
frequency bin and (4) a thin screen model with z = 4 pc with C^AL chosen to match the power at the lowest frequency bin. 
A scintillation speed uiss = 50kms"' is assumed. Spectra matched to the second lowest point are shown because they exhibit 
more power at high frequencies. However, even for this conservative choice they do not replicate the observed power at high 
frequencies. In all models all of the source flux density is assumed to undergo scintillation. For the thin-screen model a point 
source is assumed; a source of finite size further decreases the power at high frequency. In the extended medium model used here, 
the form of the power spectrum is valid for 0src > ^scat> so we choose 0src = 1-0 mas. However, the curves shown here differ little 
from ones in which the model is stretched by using O^tc = 0.1 mas instead. The parameters of the fits are summarised in Table 1. 
A minimum bias filter (Papoulis 1991, pp445-455) has been applied to the spectra in order to reduce the errors inherent in the 
spectral estimation, however the error bars shown are not mutually independent. 
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Fig. 8. (a) Two sample difference light curves, from 12 Apr 2003 (fast season) and 18 Nov 2003 (slow season). The rapid, sharp 
variations of these frequency-dependent intensity fluctuations, particularly on 12 Apr 2003, are well above the 1.4 mJy noise 
level. Note that the rapid frequency-dependent fluctuations are absent in the slow season, (b) Autocovariance functions of the 
temporal intensity fluctuations on 25 Jan 2004, 12 Apr 2004, 30 Jan 2005. Note that the amplitude of the function at zero temporal 
lag significantly exceeds the amplitude observed in the spectral autocovariance function (see Fig.|9j. The excess can be attributed 
to the additional influence of refractive scintillation on the temporal modulations. 



18 



Macquart & de Bruyn: Diffractive scintillation in J1819+3845 

22 Fe b 2004 




-5.0004 



JOO 300 300 400 500 60G 

J2 Apr 2004, 



0.0004 





CfcOiXi^ 


1 










■0.0002 



-[|.00Ol 




200 m 
30 Jan 2005 




-o:oofi4 



.m 



3Q0 -WO 



Fig. 9. Fits to the spectral decorrelation characteristics of the diffractive scintillation. The grey error regions incorporate both 
uncertainties due to the finite number of diffractive scintles incorporated in the fit and uncertainties in the mean spectrum of 
the diffractively scintillating component. The green and blue lines represent, respectively, the best fitting spectral decorrelation 
functions expected from a thin phase changing screen and from an extended scattering medium. Where only one line is visibile 
the fits are identical. The plots shown here are the averages of the autocorrelations from x 720 time slots for each observation. The 
observed decorrelation bandwidths (half-width at half maximum) are 170, 192 and 156 MHz for the three epochs respectively. 



